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SECTION  1 
INTRODUCTION 


This  report  is  concerned  with  the  analytic  investigation  of  the 
electromagnetic  portion  of  the  SGEMP  problem.  In  the  SGEMP  problem,  photo- 
electrons ejected  from  the  surface  of  an  object  are  source  currents  for 
electromagnetic  fields.  The  fields  at  the  surface  of  the  object  are  most 
interesting . 

We  investigate  SGEMP  phenomena  using  a perfectly  conducting  spheri- 
cal model  of  the  object.  The  investigation  has  three  purposes.  The  first 
purpose  is  to  present  some  directly  useful  data  of  a handbook  type  nature. 

The  second  purpose  will  be  to  get  a deeper  understanding  of  some  of  the 
physical  aspects  of  the  external  SGEMP  problem.  The  third  purpose  will  be 
to  try  and  g ..in  insight  into  SGEMP  problem  analysis  for  an  arbitrarily 
shaped  structure. 

The  time-dependent  analytic  solution  for  the  electromagnetic  fields 
in  space,  where  the  source  current  distribution  is  axisymmetric  and  the 
inner  boundary  is  a perfectly  conducting  sphere  has  been  developed  in 
Reference  1.  In  particular,  the  magnetic  field  at  the  surface  of  the  sphere, 
or  equivalently  the  skin  current,  is  kno.'n,  as  a function  of  time,  when  the 
spatial  currents  are  known  as  a function  of  time.  The  starting  point  for 
our  present  analysis  will  be  the  equation  expressing  the  dependence  of  skin 
current  on  source  currents  which  appears  in  Reference  1 . 
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To  obtain  numerical  data,  some  assumption  about  the  mathematical 
nature  of  the  source  currents  is  necessary.  Two  limiting  cases  of  the  SGEMP 
source  currents  will  be  evaluated.  The  first  case  is  that  in  which  electron 
motion  is  unaffected  by  the  fields  they  generate.  This  case  will  be  called 
the  linear  limit.  The  second  case  is  that  in  which  the  fields  generated 
by  the  system  are  so  strong  that  electrons  move  only  a small  fraction  of 
the  sphere  radius,  before  returning  to  the  sphere.  This  case  will  be  referred 
to  as  the  non-linear  limit.  Apart  from  their  simplicity,  these  two  cases 
may  prove  to  be,  when  astutely  combined,  an  important  approximation  to  the 
general  problem. 

In  order  to  make  the  transition  from  the  formal  solution  of  Refer- 
ence 1 to  an  actual  computation,  mathematics  is  required.  Section  2 will  be 
devoted  to  deriving  the  equations  necessary  for  the  computation.  Section  3 
will  characterize  the  linear  limit  and  non-linear  limit  approximations  for 
the  source  current.  Section  4 will  present  the  results  of  the  computation 
and  Section  5 will  discuss  the  results  and  their  relationship  to  the  purposes 
of  the  report.  Section  6 will  present  a summary  of  the  report  and  its 
conclusions . 
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SFCTION  2 

MATHEMATICAL  PRELIMINARIES 

This  section  will  be  concerned  only  with  the  skin  currents  on  the 
surface  of  a sphere.  (The  electric  field  at  the  surface,  or  equivalently  the 
surface  charge,  can  be  obtained  quite  simply  from  the  skin  current  and  source 
current  through  a time  integration.)  We  will  be  working  with  Equation  4-14  of 
Reference  1.  In  Reference  1 all  the  electromagnetic  fields  in  space  are  ex- 
panded in  terms  of  Legendre  polynomials.  Each  Legendre  coefficient  for  the 
fields  is  expressed  as  an  integration  over  the  complex  to  plane  and  the  spatial 
radial  coordinate  r'.  The  integrand  contains  the  Fourier  components  of  the 
Legendre  coefficients  of  the  source  currents.  Expression  4-14  of  Reference  1 
is  a particular  case  of  the  general  field  solution.  The  expression  is,  apart 
from  a multiplicative  constant,  the  Legendre  coefficient  of  the  magnetic 
field  evaluated  at  the  surface  of  the  sphere. 

In  this  section  we  will  first  put  Equation  4-14  of  Reference  1 into 
a more  usable  form.  We  will  then  obtain  some  relationships  which  enable  us 
to  transform  the  integration  over  the  complex  w domain  to  integrals  over 
time.  It  is  found  to  be  easier  to  do  the  time  integrals  than  the  complex 
plane  integrations.  The  time  integral  formulation  also  allows  a deeper  in- 
sight into  the  physics  of  the  problem.  In  the  final  parts  of  this  section 
we  will  find  a few  relations  which  simplify  the  time  integrations. 

Equation  4-14  of  Reference  1 is 
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wheie  K-  is  the  £ Legendre  coefficient  of  the  skin  current  (it  will  be 

r 9 

referred  to  as  skin  current  magnitude);  J,  and  J.  are  Fourier  integral 
th  ‘ £w 

components  of  the  £ Legendre  coefficient  of  the  r and  9 (spherical  coordin- 
ates) components  of  the  source  current  J In  and  j„  are  spherical  Hankel  and 
Bessel  functions  respectively;  R is  the  radius  of  the  sphere;  c is  the 
velocity  of  light;  ui  is  the  frequency  or  the  Fourier  transformation  and  t is 
time.  We  now  wish  to  combine  the  two  terms  in  Equation  1.  If  we  let 
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where  the  argument  of  the  spherical  Hankel  and  Bessel  functions  isua/cR,  then 
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Using  Equations  4-4  and  4-7  of  Reference  1 it  is  easy  to  show  that 
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then,  from  Equations  1,  4 and  6 we  have 

OO  00 
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Equation  7 is  the  form  of  Equation  1 we  sought.  We  now  derive  a convenient 
form  for  !f.  From  Reference  2 Equations  I - 10  to  1-12  it  is  easy  to  show  that 

■ , . , e+1  . , l 
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where  Y ( £ , m ) is  given  by  Equation  1-15  of  Reference  2: 

*(£+n.)lf:|(£2-s2) 

Y(£,mi  = — , m > 0 1 

2mm!  I 

Y(£,0)  = 1 > (9) 

0 2 2 I 

W*  - > -=  J • ) 

Having  found  a form  for  hi,  in  J\  we  consider  L.  From  8 we  find 
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By  redefining  the  summation  index  in  the  second  term  of  Equation  10, 

_ . -£ 

factoring  out  (kR)  , multiplying  the  entire  bracketed  expression  by  - i / - i 
and  using  Equation  1-45  in  Reference  2 we  find  that 
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The  expression  in  brackets  of  Liquation  11  is  a polynomial  of  degree  t,  + 1 
in  kR.  The  roots  of  this  polynomial  determine  the  electric  modal  frequencies 
of  the  sphere.  (3/3R  (Rh*)  occurs  in  the  denominator  of.'/’.)  If  we  represent 
the  q root  of  the  polynomial  of  degree  l + 1 by  then 

njlj  (kR  - pj)  * (i)1*1*  V(t.«) 


»?0  (fg:;; . «m 


Substituting  Equation  11  into  Equation  4 and  the  resultant  into 
Equation  6,  together  with  Equations  8 and  then  12  we  find  the  needed  form  for  1/’. 
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(Eventually  we  want  to  write  a dimensionless  form  for  the  skin  current  and 
Equation  16  is  a start  in  that  direction.)  We  now  wish  to  combine  Equations 
16  and  7.  Same  definitions  are  useful.  If  we  define 
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then  using  the  definition 
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in  Equation  18,  and  substituting  the  resultant  equation  together  with  Equation 
16  into  Equation  7 we  find  that 
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where  we  have  substituted  the  variable  x as  defined  in  Equation  15,  for  r’ 
in  Equations  7 and  18.  Equation  20  is  almost  the  usable  form  of  the  skin 
current  that  we  are  seeking.  One  of  the  steps  remaining  is  to  reduce  Equa- 
tion 18,  as  used  in  Equation  20,  to  an  integration  over  time.  The  following 
achieves  that  end. 

In  order  to  simplify  Equation  18  we  consider  the  following  inte- 
gral transform: 


f(t)  = 2^ f g(w)e"lwtdo>  • (21) 

-oo 

Taking  successive  derivatives  of  Equation  21  with  respect  to  t we  find  that 

CO 
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(We  assume, for  the  physical  functions  we  are  dealing  with,  that  all  the 
operations  defined  in  this  section  are  valid.)  Thus  we  see 
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Using  Equation  23  in  Equation  18,  where  we  substitute  g,  (r')  for  g (-')  (the 
time  function  f will,  as  a consequence  of  this  substitution,  be  subscripted 
with  an  £)  in  expression  23,  we  find 
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Interchanging  the  order  of  integration, while  noting  that  are  in  the  lower 
half  of  the  imaginary  portion  of  the  w plane,  we  integrate  Equation  24.  Using 
a contour  in  the  lower  half  of  the  ui  plane  for  (t  - t ’ ) > 0 (the  integral 
is  zero  for  (t  - t')  < 0)  we  find,  with  the  theory  of  residues  that 
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Equation  25  can  be  further  reduced  by  same  knowledge  of  the  roots  a.’ . The 
following  characterization  will  be  sufficient  for  our  purposes.  For  odd  l 

(an  even  degree  polynomial  determines  the  roots)  the  roots  occur  in  pairs 

i i 

so  that  if  to'  is  one  root  - to  * will  be  also.  These  roots  describe  a damped 

q q 

oscillation.  For  even  a root  can  occur  which  has  no  real  part.  This 
root  corresponds  to  damping  with  no  oscillation.  Assuming  that  to^  and  . * 
are  a pair  of  roots  we  look  at  the  form  of  the  first  two  terms  in  the  sum 
of  Equation  25: 
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1 V q=3  1 q 


anfs(t') 
3t ,n 


dt ' 


• H*, 

+ 101^  t 


„U1  V 
I!  , (-to, 
q = 3 1 


f ~ ' ^ 1 

7 J° 


<L*  , 
10),  t ’ 


dt'" 


dt  ’ 


)• 


(2b) 
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where  RP  in  Equation  26  means  the  real  part.  For  > 5,  one  ol'  the  products 
in  the  denominator  of  the  first  term  of  Equation  26  is 


nZ  , Z I l*. 

0.  - fur  - to  ) (uq  + u,  ) 

lq  1 q 1 q ’ 


(27) 


Z . 


if  the  root  ^ is  one  of  a pair  of  roots.  The  corresponding  product  in  the 


denominator  of  the  second  term  is 


lq  = (‘UT  * V('U1  + Uq  1 ' 


(28) 


It  is  clear  that 


0,  * = F, 
lq  lq 


If  Z is  odd  and  all  the  roots  occur  in  pairs  it  can  be  seen  from  Equation  28 
and  an  examination  of  Equation  26  that  the  two  terms  in  Equation  26  arc 
complex  conjugates  Of  each  other.  Under  these  circumstances  Equation  25 
becomes 

t 


A*(t)  = (i),,+“ 
n 


•>  1 


? V 1 r 

. , i t+i  rrrr  J dt  c J 

J = 1’2UjVl.2|Diq^ 


af(t'-t)  3nf  , (t ' 1 

1 k, 

sin 

3t ,n 


( 


J/<1 

Z*\ 


u.(t'  - t)  - 


Y e.  | 

1=1,2  Jq  > 


29) 


j/q 


where  we  have  defined 


Z . Z 

lx  I.  = -10  . 

J J 


J 


(50) 


U?-  = | | e+'  jq 

jq  ' jq1 


(311 


and  where  the  "2"  in  the  summation  and  multiplication  symbols  indicates 
that  only  one  root  of  a pair  is  to  be  used.  One  can  obtain  formulas  similar 
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to  Equation  29  if  all  the  roots  do  not  occur  in  pairs.  For  the  source  cur- 
rents of  this  report.  Equation  29  is  sufficient,  however.  Equations  20  and  29 
indicate  that  for  each  component  of  the  skin  current  K?  there  are  Z + 1 terms 
to  evaluate  for  each  of  the  (£  + I ) / 2 pairs  of  frequencies.  Z terms  contain 

a time  derivative  of  the  current  source  (in  Equation  29  f (t)  represents 

x*  9 0 ~ 

Jjl  + Br^"  *n  t*le  actua*  physical  system).  The  highest  order  time 

derivative  is  Z. 


By  integrating  by  parts  in  Equation  29  we  can  reduce  the  number  of 
separate  time  integrals  that  are  necessary.  In  the  computation  of  the 

following  sections  it  is  found  convenient  to  express  the  integral  of  f.  (t') 

2 ^2 

and  the  integral  of  3f.(t')/3t'  in  terms  of  the  integrals  of  3 f0(t')/3t' 

3 3 * 

anc*  ^ Ct ' ) / ’ . We  now  transform  the  time  integrals.  Defining 


«.+  l 

E 

q=l,2 


dZ 

jq 


(32) 


in  Equation  29,  we  have 


I*j 


(t) 


/dW 


o(t'-t)  3 f (t 1 ) 
3tn 


sin(a(t'-t)  - 0)  = 


t i 

9_f(tl  sin0  . f 3"  fCt')  ca(t'-t)(o  sin(a(f 

at"'1  J at11'1 


t) 


0) 


+ acos(a(t'  - t)  - 0))dt'  , (33) 

where  we  have  assumed  that  for  t < 0 the  function  and  all  its  derivatives 
are  equal  to  zero.  We  have  dropped  the  subscript  Z on  f(t)  and  the  super- 
script Z and  subscript  j on  a,  a and  0 for  convenience.  Making  a partial 
integration  in  Equation  33  we  have 
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Iij(t)  = - Sine  - (acose  - asin6) 

n 3tn-l  3tn-l 


t _ 

r 9n~2 


1 )e° ^ ^((a2  - a2)sin(a(t'  - t)  - 0) 


+ 2aacos(a(t'  - t)  - 0))dt'  . (34) 

By  combining  Equation  33  with  n = 2 and  Equation  34'  with  n = 3 in  such  a 

way  as  to  remove  the  integral  of  cos(a(t’  - t)  - 0)  we  find 

dt 

1^  = I-  - f(t)  sin0  - (asin0  + acos0) 

1 |w|2  ( 9t2  3t 

- I*J  - 2al§j } , (35) 


where 


1-1-2  2 
w = a + a . 


Similarly  by  conbining  (33),  with  n = 1,  and  Equation  34,  with  n = 2, 
together  with  Equation  35  we  find 

Tij  = -L-U 


Iwl2  !w|2  9t2 


s -) 


- F(t)[osin0  + acos0] 


2a  Ti  j 

*rfV* 


Equations  35  and  36  are  those  necessary  for  our  purposes.  As  a finale  for 
this  section  we  put  Equation  20  into  a dimensionless  form.  The  dimensionless 
form  is  convenient  for  comparing  results. 

Beginning  with  Equation  29  and  defining  a dimensionless  time  by 
t~  = t/x  , (38) 


where  x is  defined  by  Equation  14  we  find  that 


Aj(t)  = (i)n+ 


, „ .1+1  1 r I O^ft'-t 

w "*£  5V*1  ,|Dl  I /dt,'c 

J-l  .2  ) q= 1 > 2 jq1  < 
j^q 


Lsr 

3t'n 


3 lf0  Ct:’)  _ n 

sin(aj(f  - t)  - 6p)  , 


(39) 


where 


and 


Z . 


-z 

a.  = 
3 

Z 

xa  j , 

(40) 

-Z 

a.  = 
3 

Z 

xa.  , 

(41) 

D*  = 

(P^  - P^) (P^  - 

0* 

P*  ) , (42) 

39 

j a j 

q 

defined 

by  Equation  19. 

i 

It  snould  be  noted  that  0.  does  not 

1 

-Z 


-Z 


change  with  this  transformation  and  a.  and  a.  are  the  imaginary  and  real  parts 

Z ^ ^ 

of  pj  respectively.  If  we  normalize  f(t)  where  we  define 


fos 


J f(t)dt  = x J f (t)dt  = xf 


0 


(43) 


f(t)  = f(t)/fQ  , 


(44) 


then 


! 
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where 


A*(t)  = (i)n+2(x)A'nf0A^(t)  , 


1 i r ~ / a: 
S —i  £+ 1 r / ^ ' le  ^ 

j = l ,2  a ir  ; D7  ' 

J q-l ,2  jq1  0 


af(t’-t) 


_£_  _ £. 

sin(aV'  (t 1 - t)  - eh  . 

9t’  3 3 ' 


Substituting  (45)  into  (20)  where  we  make  the  correspondence 

V*')  - J£(t' ,r')  + gp-(r'Jg(t’'  ,r'))  , (47) 

— _ 
so  that  actually  A^(t)  also  contains  the  variable  r',  besides  t,  we  find  that 


K£(t)  = - cf0  Y.  (i)£+1Y(^,m)  f dx 
m=D  J 


, (48) 


where  we  have  used  the  definitions  expressed  by  Equations  14  and  15.  Equation 
48  together  with  Equation  46  are  the  dimensionless  forms  we  seek. 
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SECTION  3 


SOURCE  CURRENT  APPROXIMATION 


In  this  section  we  will  discuss  the  linear  and  non-linear  limits 
for  the  source  currents.  The  end  of  this  discussion  will  be  the  arrival  at 
a set  of  expressions,  to  be  evaluated,  which  include  all  the  parameters  of 
the  specific  source  currents  to  be  used. 

Two  approximations  are  shared  by  both  the  linear  and  non-linear 
limits.  These  approximations  are  common  to  many  "Maxwell  Solver-Particle  Pusher" 
computer  code  simulations  of  SGEMP.  As  usual,  these  approximations  are  made 
to  simplify  the  labor  involved  in  the  calculations, and  also  because  their 
affect  on  the  physics  of  a real  problem  is  not  believed  to  be  too  appreciable. 

The  first  approximation  assumes  that  electrons  are  emitted  only  perpendicular 
to  the  emitting  surface.  That  is 

J6  = 0 , (49) 

for  our  sphere  problem.  More  accurate  computer  calculations  of  SGEMP  effects 
assume  electrons  are  emitted  with  a cosine  distribution, with  respect  to  the 
normal  to  the  surface.  Although  detailed  studies  have  apparently  not  been 
performed,  the  opinion  of  some  workers  in  the  field  is  that  simulation  with 
the  two  types  of  emission  yields  similar  results.  An  upper  bound  of  30 
percent,  on  the  difference  in  magnitude  of  skin  currents  generated  by  these 
two  types  of  emission  is  sometimes,  unofficially,  quoted. 

A second  approximation  is  that  emission  from  all  points  on  the 
emitting  half  sphere  begins  simultaneously.  This  approximation  is  expressed 
as 


18 


(50) 


Jr  = H(ir/2  - 0)p(r,t)  , 

where  H is  a step  function  and  p is  a function  of  the  radial  spatial  coordin- 
ate and  time.  The  range  of  0 variable  for  the  half  sphere  that  is  emitting 
is  0 < 0 < it/2.  In  a typical  physical  problem  a plane  wave  of  X rays  strikes 
a point  of  the  sphere  (9  = 0);  all  the  points  of  the  half  sphere  begin 
emitting  after  a time  R/c  (the  time  it  takes  light  to  cross  a half  sphere). 

For  those  sources  whose  effective  time  duration  is  much  larger  than  R/c, 
Equation  50  is  clearly  a valid  approximation.  For  those  sources  who's 
effective  time  duration  is  comparable  or  smaller  than  R/c,  the  validity  of  the 
approximation  is  not  certain. 

A direct  consequence  of  Equation  50  is  that  the  electromagnetic 
fields  may  be  expressed  in  terms  of  a series  of  odd  order  Legendre  polynomials . 
The  proof  is  simple.  Expanding  H(tt/2  - 0)  in  terms  of  Legendre  polynomials 
we  have: 


H(tt/2  - 9)  = £ a£P£(cos0)  , 


L=0 


(51) 


where 
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^ / ' 


*«,  " 2 / *Vx;)dx  ’ 


(52) 


and 


x = cos9  . 


(53) 


From  the  orthogonality  properties  of  the  Legendre  polynomials  (P^(x)  is  a 
constant)  we  know  that 

0 1 


J P£(x)dx  ♦ J P£(x)dx  =0  , I?! 
-1  0 


(54) 


but,  since  P ^ (x)  is  an  even  or  odd  function  of  x when  1 is  even  or  odd. 
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U 1 

/ PJl(x)dx  = C-l/f  P£(x)dx  . 


(55) 


-1  0 

From  Equations  55  and  54  it  is  clear  that  a^  in  Equation  52  is  equal  to  zero  for 
even  (£  / 0) . Since  we  need  a^  for  our  computations  we  now  find  their  numerical 
values.  Noting  that’ 


/ 


0 

P, (x)dx 


_ (2n  - 1)!!(-1) 


n+1 


2n+1(n  + 1).' 


(56) 


where 


_ £ - 1 

n = r — 


(57) 


we  have,  after  substituting  Equation  56  into  Equation  52 


_ (4n  + 3)  (2n  - 1)  ! ■' 
n+2  fn  + 1) ! 


(-D 


n+1 


(58) 


In  this  report,  for  both  the  linear  and  non-linear  limit,  we  will 
be  concerned  only  with  the  first  two  terms  in  the  infinite  series  for  the 
skin  current,  that  is, the  terms  labeled  by  £ = 1 and  3 (n  = 0 and  1).  From 
Equation  58  we  see  that 


a1  = - 3/4  , 


a2  = + 7/16  . 


(59) 

(60) 


Looking  at  Equation  46  we  see  that  we  need  the  numerical  values  of  |D. 

£ — £ — £ —Si  — £ J 0 

0.,  a.  and  a.,  when  £ is  equal  to  3 and  a.  and  a.  when  £ is  equal  to  1,  for 
J J J J 1 

our  computations, with  either  the  linear  or  non-linear  limits  for  the  sources. 

We  begin  with  £ = 1 and  find  the  needed  values. 

The  equation  for  a1  and  o1  (the  subscript  j is  suppressed  since 
only  one  pair  of  roots  exists)  is  given  by  Equation  12  where 
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p = kR  . (61) 

Using  the  definition  of  Y(X,,m)  provided  by  Equation  9 we  see  that  the  solution 


-1+ip+p  =0, 


provides  the  roots 

a 1 = V3 77  , 

(63) 

a1  = 1/2  . 

Equation  62  and  its  solution  is  the  same  as  those  of  Reference  4.  The  equation 
describing  the  roots  for  £ = 3 is  obtained  in  the  same  way  as  the  i = 1 
equation.  It  is 

p4  + i6p3  - 21p2  - i45p  + 45  = 0 . (64) 

The  solution  of  Equation  64  is  found  by  numerically  solving  a cubic  equation. 
One  eventually  finds  that 

a3  = .870569  , a3  = 2.75785  , 

(65) 

a3  = 2.15713  , a3  = .842862  , 

— £ 

to  six  places.  Using  the  definitions  expressed  by  Equation  42,  where  D. 

■ — o | i 0 & • ^ 

is  equal  to  |Dr  |e  jq,  we  find,  after  reference  to  Equations  19,  30,  40 

3 q 

and  41,  that 

|D3iI  = |D32|  = 8.87527 


0j2  = -2.88081  radians 


021  = .955811  radians  . 


2) 


LINEAR  LIMIT 


As  inonoenergetic  electrons,  of  velocity  v,  leave  the  sphere  their 

2 

flux  decreases  by  a factor  of  (R/r1)  , where  R is  the  sphere  radius  and  r' 

2 

is  the  radial  coordinate.  If  f(t)(2TrR  ) is  the  total  number  of  electrons 
leaving  the  half  sphere  at  time  t then  that  same  number  appears  at  the  half 
sphere  located  at  r'  but  at  a time  t - (r'  - R')/v)  later5.  That  is  from 
Equations  50  and  51 

= a£(R/r')2f(t  - (r ' - R)/v)  . (67) 

Equation  67  expresses  the  linear  limit,  source  approximation.  We  wi’’  want 
to  calculate  and  plot  the  skin  currents,  due  to  a source  of  this  form,  for 
a number  of  different  velocities.  With  this  end  in  mind  we  substitute 
Equations  49  and  67  into  Equation  47  and  the  result  into  (46).  We  find, 
after  utilizing  Equations  35,  36,  37,  63,  65,  and  66  and  after  doing  some 
laborious  numerics  calculations,  that 


Ki (t)  = cfn( 


- "4>  /%  ((<- 1 * r)Fi3  - f!2)  =r 


. an . (.  ! ♦ i,  , 

x * 3t  f 


(68) 


and 


K3(t)  - - cfQ(16) 


00 

fdM 


!29424  - ^8771.  + ,^86047 


) 


»(■ 


+ .040855  - 


.073690 


♦ F*  (.  776544  - + ■ 8,75314  \ 

X / xZ  x3  ' 


+ ^22  V 245130 


.124222 


- .048510 


,)  + I MlZi  /JL  .luif, 

' 3 3t  l „2  ,.3 1 3 n"  ..3 


f(y) 


OiZl  /•  032 
^2  l ; 

3t  ' x 


313  .061551 


)]• 


(69) 
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where 


y = t - (x  - 1)(1  + y'1)  , y = v/c  , 


(70) 


Fjk(X'tJ  5 fdt 


’ f(f)e 

at1 


0j(t’*y)  . 


sin (ar(t'.  y)  - 07)  , (71) 


and 


E0  • /' 


f„  = I f(t)dt  , 
0 


(72) 


(the  total  amount  of  charge/cm2  leaving  the  conducting  surface)  where 

f(t)  5 (f(t)T)/fQ  . (73) 

The  in  expression  69  has  six  place  accuracy.  Note  that  f(t)  is  defined  by 
Equation  67  and  also  that  the  definitions  expressed  by  Equations  IS  and  38  have 
been  used  in  deriving  Equations  68  and  69.  The  other  symbols  used  under  the  in- 
tegral sign  of  Equation  71,  to  be  used  in  Equations  68  and  69,  can  be  obtained 
from  Equations  63,  65  and  66,  where  we  note,  however,  that  9*  is  equal  to  zero. 
Equation  70  defines  a dimensionless  retarded  time  y with  an  effective  velocity 
of  retardation  ce,  where 


Cg  _ 1 

— = 1 + Y • (74) 

Both  the  time  for  a charged  particle  to  get  from  R to  r’,  and  the  time  for 

light  to  travel  from  r'  to  R are  included  in  y.  The  reason  y occurs  in 

expressions  (68),  (69)  and  (71)  rather  than  t - (x  - 1)  can  be  seen  by 

f r 1 -R1 

noting  that  the  Fourier  transform  of  f(t  - given  by 


-OO 


f(t)dt  , 


or  is  equal  to  the  Fourier  transform  of  f(t)  multiplied  by  the  factor 
e+iw((r  R)/v).  This  factor  has  the  same  form  as  the  exponential  factor 
eik(r'-R)  -n  Equatjon  antj  so  (r'-R)/v  would  appear  in  the  subsequent 
equations  of  Section  2 in  the  same  way  that  (r'-R)/c  appears. 
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More  useful  forms  of  Equations  68  and  69 — from  a computational 

£ 

standpoint — can  be  found  by  partially  integrating  the  F.^  terms  in  these 
equations  while  also  expanding  the  sine  and  cosine  terms  in  Equation  69 
The  resulting  equations  are 


Kj(t)  = cf 


3 fdx  ( f(> 

0 4/  2 ' x 

J x 


+ n . I)  M (y) 


3t 


- v? 11  * c|2  * a-i)sj,i  , 


xJ  "12» 


(68') 


and 


■ - <*„&>  ff2  (l  w 7 * I ^7  - 7) 


451492 


.667222 


„3  l .130158  . 

+ I-.  029531  2 + 


,^P),C32( 

. 068899^  „3  / 

„3  ) Q22\- 


237098 


.501517  .307195 


237098 


+ -307195  V 

^ / 


.167983  .0261381 


(6  9’) 


where 


w 


*(111 


y)  _£  _ 

sina. (t 1 -y)di ' , 


0 3t’k 


and 


y k _ ^(t’-y) 

_ r 9 f(t’) 


/ 


9t,k 


cosol  (t  '-y)dt ' . 


(71’) 


In  general  photoelectrons  leaving  the  spherical  surface  are  not 
monoenergetic  but  are  described  by  a spectrum  of  velocities.  If  f(t)  is 
known  for  each  of  the  different  electron  velocities  then  K^Ct)  can  be  summed 
for  each  of  the  individual  velocities  to  obtain  the  effect  of  the  total 
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spectrum.  If  all  the  photoelectrons  have  the  same  time  history,  as  it  is 
usually  assumed  in  SGEMP  work,  the  total  skin  current  calculation  is  simpler. 

NON-LINEAR  LIMIT 


In  this  limit  we  assume  the  fields  generated  are  so  strong  that 
electrons  ejected  from  the  surface  move  only  a small  fraction  of  the  radius 
before  returning.  Under  these  circumstances  the  spatial  integration  in 
Equation  48,  the  integration  over  x,  is  basically  that  of  a delta  function 
at  x = 1,  that  is  the  integrand  is  evaluated  at  x equal  to  1.  The  physical 
quantity  that  results  from  the  spatial  integration  can  be  viewed  as  the  time 
rate  of  change  of  a dipole  moment  per  unit  area.  The  dipole  moment  has  a time 
history  f(t).  As  in  the  linear  limit.one  can  view  the  total  skin  current 
as  coming  from  a sum  of  the  effects  of  different  dipole  moments  produced  by 
several  groups  of  monoenergetic  electrons  with  different  time  histories. 

One  could  also  view  the  dipole  moment  of  the  spatial  integrations  as  the 
total  dipole  moment  of  several  groups  of  electrons  with  the  same  time  history. 
This  latter  view  is  closer  to  the  usual  SGEMP  assumption.  Integrating  an  x 
in  Equations  68'  and  69'  as  though-  tire  integrand  contained  S(x  - 1)  we  find 
the  skin  current  is  f(t)  (see  Equation  73  for  a definition  of  f(t)): 


Together  with  Equations  70,  71'  (setting  x = 1),  63  and  65,  Equations  75 
and  76  are  those  needed  for  the  non-linear  limit.  We  will  use  these 
equations  to  compute  skin  currents  for  various  f(t). 

it  should  be  noted  that  in  arriving  at  Equations  75  and  76  we 
assumed  that  the  source  currents  had  the  form  S(x  - 1).  This  assumption 
becomes  invalid  for  1 greater  than  some  finite  integer  "N,”  because  the 
wavelengths  of  the  oscillations  will  become  comparable  to  and  smaller  than 
the  dipole  thickness.  To  make  the  series  of  convergent, in  a mathematical 
sense,  for  the  non-linear  limit, an  integration  over  x would  have  to  be 
performed  for  l > N. 
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SECTION  4 
RESULTS 


We  begin  first  with  the  non-linear  limit  since  it  is  simpler  and 
consequently  demonstrates  the  implications  of  the  theory  more  clearly.  To 
do  the  computations  it  remains  only  to  choose  a time  history  (pulse  shape) 
and  substitute  it  into  Equations  75  and  76.  A simple  time  history,  although 
not  necessarily  adequate,  physically,  is  a triangular  pulse  which,  starting 
from  t = 0,  reaches  its  maximum  in  t^/2  and  is  over  at  t^  (t^  will  be  referred 
to  as  pulse  length).  In  terms  of  the  dimensionless  time  and  a normalized  pulse 
we  have 

f(t)  = 6t  0 < t < tf/2 

f(t)  = B(tf  - t)  tf/2  < t < tf  (77) 

f(t)  =0  t f < t < 0 , 



= 3(H(0)  - 2 H(t  - t f/2)  + H(t  - t~))  , (78) 

i r 


— = 6(6(0)  - 26 (t  - tf/2)  + 6 (t  - t ))  , (79) 

3t^ 

^4  = 6 4 6(0)  - 2 4 6(t  - tf/2)  + 4 6(t  - t,))  . (80) 

3t3  3t  3t  * 3t  f 

Where 

4 = tf / x , (81) 

6 = (tf/2)'2  , (82) 
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H is  a step  function  and  6 is  a delta  function.  Substituting  Equations  77 
through  80  into  Equations  75  and  76  and  plotting  the  results  for  various 
values  of  tf  we  obtain  Figures  1,  2 and  3.  We  have  divided  K,  by  the  factor 
c fg  and  also  by  the  factor  gt^.  to  obtain  the  ordinate  values.  Figure  1 
is  a plot  of  Equation  75  for  the  above  defined  f(t).  As  will  be  explained 
in  the  next  section  under  "Low  and  High  Frequency  Limits,"  the  term  c ~ f(^f(t), 
in  Equation  75,  is  the  quasi-static  limit  of  the  £ = 1 component  of  the  skin 
current.  In  Figure  1 the  quasi-static  limit  appears,  for  any  pulse  length, 
as  a triangle  of  height  .5  with  a width  of  1.0.  The  plots  demonstrate  how 
the  actual  skin  current  differs  from  the  quasi-static  limit  for  the  t values 

l * 

of  .1,  .5,  1,  5 and  10.  Equation  71'  defines  the  term  S which  appears  in 

— £ 

Equation  75.  Setting  ck  equal  to  zero  in  that  expression  yields  the  graph 
labeled  (1,0)  in  Figure  1— the  pulse  length  t^  is  1 and  there  is  no  damping. 

T.. . plots  of  Figure  2 are  similar  to  those  of  Figure  1.  Two  curves  are  plotted 
along  with  the  quasi-static  limit.  Both  are  plots  of  Equation  75  with  a 
pulse  length  of  12.  One  plot  contains  the  damping  in  the  s|_,  term  the  other, 
labeled  (12,0)  does  not. 

The  quasi-static  limit  for  Equation  76  is  the  term  -c  f.  f(t). 

Id  U 3 

Using  the  triangular  pulse  function  defined  above  for  f(t)in  Equation  76 
we  find  values  for  the  plots  of  Figure  3.  In  Figure  3 the  ordinate  values  of 

skin  current  magnitude  are  divided  by  -c  fggtj.  while  the  abscissa  is  t/F^.. 

As  in  Figures  1 and  2,  the  quasi-static  limit  for  £ = 3 appears  as  a tri- 
angle of  height  .5  and  width  1.  The  pulse  widths  considered  were  t^  = .1, 

. 5 , 1 , 5 and  1 0 . 

The  triangular  pulse  defined  by  Equations  77  is  discontinuous. 

Figure  4 is  a plot  of  the  £ = 1 skin  current  response  for  a smooth  function 
with  a continuous  first  derivative  at  the  beginning  and  at  the  end  of  the 
pulse.  The  time  history  is  defined  as  follows: 
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Figure  2.  Z = 1,  non-linear  limit,  triangle  excitation.  Normalized 
skin  current  magnitude  vs.  normalized  time. 
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lear  limit,  double  triangle  excitation, 
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angle  excitation. 
:ed  time. 


Two  pulse  lengths,  t£  = 1,  12,  are  compared  with  the  quasi-static  limit  in 
Figure  4.  It  is  of  interest  to  compare  the  deviations  from  the  quasi-static 
limit,  for  the  corresponding  pulse  lengths  in  Figures  1,  2 and  4.  Another 
type  of  pulse  shape  can  be  defined  by  a double  triangle  as  follows: 

f(t)  = Bt  0 < t < t{/2  , 

f(t)  = B(tf  - t)  tf/2  <t  <ytf  , 

_ _ (85) 

f(t)  = B(2tf  - t)  t£/2  < t < 2tf  , 

f(t)  =0.0  2tf  < t , 

where 

6 5 (tf/2)2  , (86) 

and  fg  in  Equations  75  and  76  is  the  time  integral  of  the  actual,  non- 
normalized,  time  function  over  the  range  0 < t < xt£/2.  Plots  of  skin  cur- 
rent magnitudes  produced  by  the  latter  time  history  appear  in  Figures  5 and 
6.  They  are  plotted  in  a manner  similar  to  the  plotting  of  the  single  tri- 
angle time  history. 

Plots  of  ft  = 1 and  ft  = 3 skin  current  magnitudes,  for  the  linear 
limit,  appear  in  Figures  7 through  28.  The  skin  current  values  for  Figures 
7 through  26  were  calculated  by  substituting  the  triangular  time  history, 
defined  by  Equations  77,  into  Equations  68' and  69' and  performing  a spatial 
integration  in  x.  Five  pulse  lengths,  t£,  and  five  v/c  ratios  are  used. 

The  pulse  lengths  are  .1,  .5,  1.0,  5 and  10.  The  v/c  values  are  .063,  .136, 
.186,  .239  and  1.0.  The  ordinate  values  for  the  ft  = 1 plots  are  the  skin 

3 

current  divided  by  c j f.;  the  ordinate  values  for  the  ft  = 3 plots  are  the 

4 7 4 1 

skin  current  divided  by  -c  -jy-  fg  j j.  Sometimes  the  ordinate  values  of 
skin  current  magnitudes  are  also  multiplied  by  10  or  102.  As  in  Figures  1 
through  6 a comparison  is  made  wirh  the  quasi-static  limit.  The  skin 
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Figure  7.  I = 1 , linear  l.imit,  triangular  pulse  excitation  (tf  = .1). 
Skin  current  magnitude  vs.  normalized  time. 

1 = quasi-static  limit,  v/c  = .186;  1 = complete  solution,  v/c  = .186 

2 = quasi-static  limit,  v/c  = .136;  2 = complete  solution,  v/c  = .136 

3 = quasi-static  limit,  v/c  = .063;  3 = complete  solution,  v/c  = .063 
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28.  l = 3,  linear  limit,  triangular  pulse  excitation  (tf  = .1). 

Skin  current  magnitude  vs.  normalized  time. 

” quasi-static  limit,  v/c  = .186;  1 = complete  solution,  v/c  = .186 

— quasi-static  limit,  v/c  = .136;  2 = complete  solution,  v/c  = .136 

quasi-static  limit,  v/c  = .063;  3 = complete  solution,  v/c  = .063 
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Figure  9.  I = 1 , linear  limit,  triangular  pulse  excitation  (t1 
Skin  current  magnitude  vs.  normalized  time. 

T = quasi-static  limit,  v/c  = 1 ; 1 = complete  solution,  v/c  = 

2 = quasi-static  limit,  v/c  = .239;  2 = complete  solution,  v/c 


Figure  11.  SL  = 1 , linear  limit,  triangular  pulse  excitation  (t 
Skin  current  magnitude  vs.  normalized  time. 

1 = quasi-static  limit,  v/c  = .186;  1 = complete  solution,  v/c  : 

2 = quasi-static  limit,  v/c  = .136;  2 = complete  solution,  v/c  ; 

3 = quasi-static  limit,  v/c  = .063;  3 = complete  solution,  v/c  = 


Figure  12.  i = 3,  linear  limit,  triangular  pulse  excitation  (tf  = 
Skin  current  magnitude  vs.  normalized  time. 

— - --quasi-static  limit,  v/c  = .186;  1 = complete  solution,  v/c 

quasi-static  limit,  v/c  = .136;  2 = complete  solution,  v/c 

quasi-static  limit,  v/c  = .063;  3 = complete  solution,  v/c 
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Figure  14.  1 = 3,  linear  limit,  triangular  pulse  excitation  (t^  = .5). 

Skin  current  magnitude  vs.  normalized  time. 

T = quasi-static  limit,  v/c  = 1;  1 = complete  solution,  v/c  = 1 
2 = quasi-static  limit,  v/c  = .239;  2 = complete  solution,  v/c  = .239 
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Figure  15.  I = 1,  linear  limit,  triangular  pulse  excitation  (tf  = 1.0). 
Skin  current  magnitude  vs.  normalized  time. 

1 = quasi-static  limit,  v/c  = .186;  1 = complete  solution,  v/c  = .186 

2 = quasi-static  limit,  v/c  = .136;  2 = complete  solution,  v/c  = .136 

3 = quasi-static  limit,  v/c  = .063;  3 = complete  solution,  v/c  = .063 
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00 


= 3,  linear  limit,  triangular  pulse  excitation  (t^  = 1 
kin  current  magnitude  vs.  normalized  time. 


- quasi-static  limit,  v/c  = .186;  1 = complete  solution,  v/c 

— quasi-static  limit,  v/c  = .136;  2 - complete  solution,  v/c 

■■  quasi-static  limit,  v/c  = .063;  3 - complete  solution,  v/c 
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Figure  18.  £ = 3,  linear  limit,  triangular  pulse  excitjtion  (t*  = 1.0) 

Skin  current  magnitude  vs.  normalized  time. 

1 = quasi-static  limit,  v/c  = 1;  1 = complete  solution,  v/c  = 1 

2 = quasi-static  limit,  v/c  = .239;  2 = complete  solution,  v/c  = .239 


Figure  19.  I = 1 , linear  limit,  triangular  pulse  excitation  (t^  = 5. 
Skin  current  magnitude  vs.  normalized  time. 

■■  ■ - - — quasi-static  limit,  v/c  = .186;  1 = complete  solution,  v/c 

quasi-static  limit,  v/c  = .136;  2 = complete  solution,  v/c 

quasi-static  limit,  v/c  = .063;  3 = complete  solution,  v/c 


Figure  20.  1 = 3,  linear  limit,  triangular  pulse  excitation  (t 

Skin  current  magnitude  vs.  normalized  time. 

f = quasi-static  limit,  v/c  = .186;  1 - complete  solution,  v/c 

2 = quasi-static  limit,  v/c  = .136;  2 = complete  solution,  v/c 

3 = quasi-static  limit,  v/c  = .063;  3 - complete  solution,  v/c 
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Figure  22.  1 = 3,  linear  limit,  triangular  pulse  excitation  (t 

Skin  current  magnitude  vs.  normalized  time. 

1 = quasi-static  limit,  v/c  = 1;  1 = complete  solution,  v/c  = 1 

2 = quasi-static  limit,  v/c  = .239;  2 = complete  solution,  v/c 


Figure  23.  si  = 1 , linear  limit,  triangular  pulse  excitation  (t, 
Skin  current  magnitude  vs.  normalized  time. 

1 = quasi-static  limit,  v/c  = .186;  1 = complete  solution,  v/c  : 

2 = quasi-static  limit,  v/c  = .136;  2 = complete  solution,  v/c  = 

3 = quasi-static  limit,  v/c  = .063;  3 = complete  solution,  v/c  = 


Figure  25.  Z = 1 , linear  limit,  triangular  pulse  excitation  (t. 
Skin  current  magnitude  vs.  normalized  time. 

1 = quasi-static  limit,  v/c  = 1;  1 = complete  solution,  v/c  = 1 

2 = quasi-static  limit,  v/c  = .239;  2 = complete  solution,  v/c  ; 


Figure  26.  1 = 3,  linear  limit,  triangular  pulse  excitation  (t 

Skin  current  magnitude  vs.  normalized  time. 

1 = quasi-static  limit,  v/c  = 1 ; 1 = complete  solution,  v/c  = 1 

2 = quasi-static  limit,  v/c  = .239;  2 = complete  solution,  v/c  = 


Statamp/cm  x (3/4  cfQ) 
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Figure  27.  £ = 1 , linear  limit,  triangular  pulse  excitation  (tf  = 1.0). 

Skin  current  magnitude  vs.  normalized  time.  No  damping. 

1 = quasi-static  limit,  v/c  = .186;  1 = complete  solution,  v/c  = .186 

2 = quasi-static  limit,  v/c  = .136;  2 = complete  solution,  v/c  = .136 

3 = quasi-static  limit,  v/c  = .063;  3 = complete  solution,  v/c  = .063 
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Figure  28.  £ = 1,  linear  limit,  triangular  pulse  excitation  (tf  = 1.0). 

Skin  current  magnitude  vs.  normalized  time. 


1 = quasi-static  limit,  v/c 

2 = quasi-static  limit,  v/c 


1;  1 = complete  solution,  v/c  = 1 
.239;  2 = complete  solution,  v/c  = .239 
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current  values  for  Figures  27  and  28  were  made  in  the  same  way  as  those  of 

_ £ 

Figures  7 through  26  except  that  ck  was  set  equal  to  zero  in  Equation  71 ' . 
These  latter  two  figures  are  made  for  t^  = 1 only.  Table  1 is  a summary  of 
two  aspects  of  the  total  information  contained  in  Figures  7 through  26.  The 

3 

maximum  values  of  the  skin  current  magnitudes  (divided  by  — cf^)  are  given 
along  with  the  rise  time  (in  units  of  t/x)  to  attain  this  maximum  value.  The 
corresponding  quasi-static  values  are  given  for  comparison. 

An  example  of  how  one  might  use  Table  1 is  the  following.  Suppose 

10  2 -8 
2 x 10  electrons/cm  leave  a sphere  of  radius  3 meters  in  10  x 10  seconds. 

We  have,  noting  that  3 x 10^  cm/sec  is  the  velocity  of  light. 


T = 


300 


3 x 10 


10 


= 10  sec 


rf 


10  x 10 


-8 


10 


-8 


= 10  . 


If  we  are  convinced  that  most  of  the  escaping  electrons  have  a v/c  value  of 
.136  then  looking  at  the  table  for  a pulse  length  of  10  we  find  under  "complete 
solution"  a rise  time  of  7.25  and  a maximum  normalized  skin  current  of  5.73  x 


10 


-2 


The  actual  rise  time  then  is 


x (7 . 25)  = 7.25  x 10~8  sec  . 
The  actual  skin  current  is 


| cfQ(5.73xlO*2)  = (.75) (3xl010) (2xl010x4.8xl0'10) (5.73xl0‘2) 

= 1.25  x 1010  *1*3*2. 

cm 

4.8  x 10  is  the  charge  of  an  electron.  (To  convert  to  amps/m  we  multiply 

_7 

by  .3  x 10  . For  most  all  problems  in  the  linear  limit  which  require  only 

the  rise  time  and  the  maximum  skin  current,  the  l = 1 solution  is  a good 
approximation.  The  table  is  convenient  for  those  problems  which  require  a 
linear  interpolation  between  table  values. 
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SECTION  5 
DISCUSSION 


In  this  section  the  results  presented  in  the  previous  section  will 
be  discussed.  To  aid  in  understanding  the  results,  reference  will  continually 
be  made  to  the  quasi-static  limit.  Apart  from  understanding  the  results, 
a point  of  view  which  we  want  to  examine  is  to  what  extent  the  quasi-static 
limit  can  be  used  as  an  approximation  to  the  actual  skin  current.  Put  in 
another  way,  we  would  like  to  judge  how  much  of  the  sphere  response  (the 
number  of  modes)  need  we  treat  in  a fully  time-dependent  manner  and  how  much 
of  the  response  can  we  legitimately  approximate  quasi -statical ly.  For  sys- 
tems which  are  not  as  highly  damped  as  the  conducting  sphere, we  would  also 
like  to  gain  some  insight  into  the  validity  of  this  viewpoint. 

The  quasi-static  or  low-frequency  limit,  for  the  skin  current  on  a 
sphere,  was  derived  in  Reference  1 by  taking  the  limiting  form  of  the  Green's 
Function  (the  same  effect  could  be  achieved  by  letting  k -»  0 in  Equation  13). 
In  our  discussion  of  low-  and  high-frequency  limits  at  the  end  of  this  sec- 
tion, we  take  the  limits  of  the  skin  currents,  which  were  obtained  by  using 
the  exact  Green's  function.  Looking  at  Figures  1 and  3 it  appears  that, 
for  t^  > 10,  in  the  non-linear  limit,  the  quasi-static  limit  will  be  a very 
reasonable  approximation  to  the  actual  skin  current  response  of  the  sphere. 
(For  definiteness  we  will  assume  that  the  quasi-static  solution  is  a good 
approximation  when,  foT  a given  "it,''  three  conditions  are  satisfied:  (1) 

the  main  peak  of  the  quasi-static  solution  differs  by  less  than  15  percent 
from  the  main  peak  of  the  fully  time  dependent  solution,  (2)  the  rise  time 
to  the  peak  of  the  quasi-static  solution  differs  from  the  corresponding  rise 
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time  of  the  fully  time-dependent  solution  by  less  than  15  percent  and  (3) 

the  secondary  peaks  of  the  main  solution  are  less  than  15  percent  of  the 

peak  of  the  quasi-static  solution.)  By  examining  Equations  68',  69' and  71' 

it  is  clear  that  the  magnitude  of  the  oscillatory  responses  (the  terms  that 

basically  cause  the  time  dependent  solution  to  differ  from  the  quasi-static 

—i 

limit)  depends  upon  the  damping  constant  a..  Figure  2 demonstrates  what 

—Si  -1 

happens  if  the  damping  is  zero  (a.  = 0)  in  the  integral  of  ciuation  71',  for 


£ = 1 and  tf  = 12. 


The  quasi-static  solution  would  no  longer  meet  our 
criteria  for  a good  approximation.  This  latter  argument  should  be  taken  as 
more  intuitive  than  exact  since  the  factors  that  multiply  the  oscillatory 

—l 

response  terms  depend  upon  the  magnitude  of  a..  Nonetheless  it  is  probably 
correct  to  say  that  strong  damping  extends  the  validity  of  the  quasi-static 
treatment.  This  is  a particularly  useful  criteria  for  picking  important 
modes  of  arbitrary  structures.  In  determining  whether  or  not  a mode  with 
period  T will  be  excited  by  a pulse  whose  length  is  t^  one  usually  judges 
that  if 


T/tf  ~ 1 


(87) 


the  mode  will.be  strongly  excited.  For  a system  having  a particular  mode  which 
damps  with  an  e folding  time  of  t^  the  mode  will  not  be  as  strongly  excited  if 

tf/fp  > n , (88) 

(n  is  a number  characteristic  of  the  mode),  in  spite  of  the  satisfaction  of 
Equation  87.  In  the  sphere  problem,  for  1 < t^  < 10,  treating  the  £ = 1 
mode  in  a fully  time-dependent  manner  while  treating  the  rest  of  the  modes 
quasi-statical ly  would  result  in  a reasonable  approximation  for  the  skin 
current.  For  t^  < 1,  in  Figure  1 the  magnitude  of  the  fully  time-dependent 
response  differs  quite  cc.’.oiderably  in  magnitude  and  time  dependence  from 
the  quasi-static  response.  As  t^  approaches  zero, both  the  £ = 1 and  £ = 3 
(Figure  3)  time-dependent  responses  approach  zero.  Investigation  of  modal 
response,  for  this  example,  in  the  pulse  length  range  tf  < 1 would  require 
examination  of  higher  modes. 
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The  oscillatory  portion  of  Equation  75' is  contained  in  the  term 
which  is  defined  in  Equation  71'.  Since  involves  a second  derivative, 
oscillations  are  stimulated  when  the  slope  of  f(t)  changes.  The  slope  of 
the  triangular  time  function  defined  by  Equations  77  to  80  changes, 
discontinuously , at  three  points.  Specifically  for  t < t^/2.  Equation  75 
yields 

Kx  = c | fQB(t  + ~ e'-5tsin  VT/2t)  , (89) 

where  use  has  been  made  of  Equations  63.  When  the  slope  of  f(t)  changes 
continuously  cannot  be  represented  as  simply  as  Equation  89.  Figure  4 
shows  the  effect  on  the  skin  current  magnitudes  when  a smooth  continuously 
changing  f(t)  is  used.  A comparison  of  Figures  4,  2 and  1 show  that  the 
effects  on  the  1 = 1 skin  currents  are  much  the  same  (relative  to  tie  quasi- 
static solution),  with  a gi^en  pulse  length,  for  either  the  smooth  or  the 
discontinuous  pulse.  A number  of  things  can  be  seen  by  inspecting  the 
simple  form  of  Equation  89.  The  oscillatory  term  begins  from  zero  at  t = 0. 

As  the  sine  term  gets  larger  the  amplitude  drops  due  to  the  damping  exponential. 
The  sine  term  reaches  its  maximum  at  T/4  (T  is  the  normalized  period  of  oscil- 
lation for  H = 1)  which  is  tt/VT  for  1 = 1 . At  t = T/4  in  Equation  89  the 
ratio  of  the  oscillatory  term  to  the  quasi-static  term  (t)  is  .25.  If  the 
damping  term  were  not  present  the  ratio  of  the  two  terms  would  be  .64.  This 
reasoning  leads  us  to  a more  precise  statement  of  Equations  87  and  88  for  the 
Si  = 1 mode  with  a triangular  pulse.  The  term  e ,5t  sin  VT/2  t reaches  its 
maximum  of  .47  at  t = 1.2, so  that  in  order  for  the  oscillatory  term  to  be  only 
15  percent,  say,  of  the  maximum  of  the  quasi-static  term  the  expression 

•15  V2  (-47)  ' (90) 


must  be  satisfied.  tf,  then  is  determined  by 


7.2  . 


(91) 


Since  tp  = 2,  we  have,  in  this  case 


j 


! 

! 
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That  is  ri  in  Equation  88  is  roughly  4.  Since  the  normalized  period  of  oscil- 
lation is  4tt/VT,  we  also  have 


T/tf  = 1.0  . (93) 

The  characteristics  of  Figures  5 and  6 are  much  like  those  of 
Figures  1 through  3 and  so  we  proceed  to  a discussion  of  the  linear  limit  and 
Figures  7 through  28.  Table  1 summarizes  some  of  the  pertinent  data  contained 
in  these  graphs.  The  table  shows  the  maximum  normalized  skin  current  rise 
time  for  the  H = 1,  complete  time-dependent  and  quasi-static  solutions, 
the  numbers  in  the  upper  left  hand  corners  of  the  table  entries,  appropri- 
ate to  the  quasi-static  solution,  are  the  percent  difference  between  the 
complete  solution  and  the  quasi-static  solution.  For  a pulse  length  t^  = 10 
both  the  rise  time  and  the  maximum  are  within  15  percent  of  each  other  for 
both  solutions.  When  tf  = 5 the  rise  times  begin  to  differ  by  more  than  15 
percent  when  v/c  > .186.  A look  at  Figures  22  and  24  shows  that  the  X,  = 3 
quasi-static  solution  is  a good  approximation  even  when  v/c  > .186  at  t^  = 5. 
If  tj.  = 1 the  difference  in  rise  times  between  the  two  solutions  in  Table  1 
is  roughly  60  percent.  T'  -■>  magnitudes,  however,  are  surprisingly  close  for 
v/c  < .239. 

Figures  27  and  28  are  plots  of  skin  current  magnitude  for  t^  = 1 
where  the  damping  has  been  set  equal  to  zero  in  the  integrals  defined  by 
Equation  71'.  These  curves  are  suggestive  of  the  differences  that  would  oc- 
cur between  the  complete  time -dependent  solution  and  the  quasi-static 
solution  for  a system  whose  modes  aren't  damped  as  strongly  as  those  of  the 
sphere.  Large  differences  occur  in  the  first  peak  and  not  just  after  the 
driving  pulse  is  over. 


a table,  such  as  Table  1,  for  both  the  £ = 1 and  £ = 3 modal  responses 
would  define  the  ranges  of  validity  of  the  quasi-static  solution.  Practically 
speaking,  it  is  doubtful  that  many  problems,  involving  a conducting  sphere, 
in  the  linear  limit,  would  require  more  than  one  or  two  fully  time-dependent 
modes  for  an  accurate  description,  over  most  of  the  interesting  time  range. 


LOW  AND  HIGH  FREQUENCY  LIMITS 


In  this  portion  of  Section  4,  starting  with  Equations  68’ and  69' 
we  will  first  derive  the  low-frequency  limit.  The  low-frequency  or  quasi- 
static limit  is  mentioned  throughout  this  section.  Using  Equations  68' and 
69’  again  we  outline  the  derivation  the  high-frequency  limit  in  the  final 
part  of  this  section.  The  high-frequency  limit  is  important  for  appreciating 
results  in  the  ncn-linear  limit  for  tf  < 1.  If  we  assume  that 

Yf » '• 

in  Equation  71  then  Equation  68’  becomes 


( 1 


U 3f(y) 

x 3t  , 


(94) 


(95) 


Analogously  Equation  69’  becomes 


,.  r 7 fdx  /f(y)  . K 3f(y)j 

K2(t)  = - cf0  48  J T lx  + (1  ' ^t  / 


(96) 


X V 'X 

If  t » (r  - R)c~*  (the  velocity  of  light  is  assumed  to  approach  infinity) 
we  also  have  t"  » (x  - 1).  If  we  also  assume  that  v/c  <<  1 then  we  can 
expand  f(y)  about  tf  - (x  - l)y  using  the  definition  expressed  by  Equation  70. 
The  expansion  is 
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f(y)  = f(t  - (x  - l)y_1)  - (x  - 1)  M.  (t  - (x  - 1)y_1 

at 

+ *?((x  - l)2)  . 

Substituting  Equation  97  into  Equations  95  and  96  we  have 

K^t)  = cfQ(J)  r-f  fCt  - (x-  1)Y_1)  7 , 

j X 

K2(t)  = - cf  (^r)  f- jf(t  - (x  - Dy'1)  ± 

J X X 


) 

(97) 


(98) 

(99) 


Equations  98  and  99  are  expressions  of  the  low-frequency  linear  limit.  By 
neglecting  the  velocity  and  assuming  that  the  integrands  in  this  latter 
expression  are  multiplied  by  6(x  - 1)  we  get  the  non-linear  limit: 


Kj(t)  =cf0(|)f(t)  , (100) 

K2(t)  = - cf0(^)f(t).  (101) 


The  high-frequency  limit  is  found  by  assuming  that 


o-(t'-t) 


= 1 + oJ(t'  - t)  , 


and  that 


cos(a^(t'-t))  = 1 , 
sin  (ctj  (t ' -t))  = a^(t'-t)  , 


(102) 


(105) 


in  Equations  68'  and  69'.  The  resulting  equations  are 

K1  (t)  = cf0(f)  f f)(*  J + Hf(t)  , (104) 

J x 

and 
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k2">  s - cf0 '??)/  fJT-"  - b ■ 

^ X X 

The  integrand  in  expressions  (104)  and  (105)  is  zero  at  x = 
in  Equation  105  have  five  place  accuracy. 
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SECTION  6 

SUMMARY  AND  CONCLUSIONS 

This  report  can,  perhaps,  be  divided  into  four  areas  of  endeavor. 
These  areas  include  the  systematizing  of  some  mathematical  operations,  an 
attempt  at  ellucidating  a physical  process,  the  presentation  of  some  useful 
data  and  the  expression  of  opinions  as  to  the  general  applicability  of  modal 
analysis  in  SGEMP. 

In  order  to  investigate  the  effect  on  surface  currents  of  a variety 
of  external,  axially  symmetric  source  conditions,  a good  deal  of  mathematical 
manipulations  are  necessary.  The  manipulations  are  required  to  make  the 
transition  from  a formal  solution  to  a usable  formalism.  The  solution  for  the 
currents  on  the  surface  of  a conducting  sphere  is  developed  for  odd  £,  where 
£ designates  the  order  of  the  Legendre  polynomial  expansion  functions.  Al- 
though not  explicitly  given  in  the  report,  the  solution  for  even  £ proceeds 
along  the  same  lines. 

Data  is  presented,  for  the  £ = 1 and  £ = 3 modes,  with  two  assumed 
types  of  spatial  source  current  configurations  and  many  different  source 
current  time  histories.  This  data  is  presented  in  graphs  and  a table  so  that, 
under  certain  restrictions,  some  estimates  can  be  made  in  practical  situations. 
The  restrictions  pertain  to  the  model  of  the  source  currents  and  how  well  an 
object,  in  practice,  conforms  to  the  geometry  of  a conducting  sphere. 

The  data  for  the  time-dependent  £ = 1 and  £ = 3 solutions  is 
compared  with  the  low-frequency  or  quasi-static  solution.  Comparisons  are 
made  because,  in  general,  quasi-static  solutions  are  simple  to  obtain  and 
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below  a certain  excitation  frequency  range  they  are  a very  good  approximation 
to  the  whole  true  solution.  The  true  time-dependent  solution  for  these  types  of 
source  current  configurations  — at  least-is  found  to  be,  basically,  the  sum  of  the 
quasi-static  solution  plus  oscillatory  terms.  The  oscillatory  terms  are 
time  integrals  containing  the  damping  and  oscillatory  frequencies  of  the 
spherical  modes.  Because  of  the  high  damping  of  the  spherical  modes,  the 
oscillatory  terms  contribute  less  than  might  be  expected  so  that  the  quasi- 
static solution  is  a valid  estimate  for  a large  range  of  source  time  histories. 
The  treatment  of  the  sphere  then,  in  most  situations  encountered  in  practice, 
would  probably  not  warrant  more  than  the  utilization  of  one  or  possibly  two 
modes  in  a fully  time-dependent  manner.  The  other  modes  > 3,  say)  could 
be  treated  quasi-statically , if  it  was  deemed  necessary  to  use  them. 

The  situation  with  the  sphere  suggests  that  the  fields  at  the  sur- 
face of  more  complicated  objects  could  possibly  be  treated  by  a hybrid 
analysis  of  a few  complete  time-dependent  modal  contributions  plus  a quasi- 
static' remainder.  Criteria  for  choosing  the  modes  to  treat  in  a completely 
time-dependent  manner,  would  depend  on  damping  rate  as  well  as  their 
frequency  of  oscillation.  Generally  speaking,  highly  damped  modes  would 
contribute  to  the  quasi-static  remainder  only.  The  integrals  that  describe 
the  oscillatory  contribution  of  a mode  are  probably  similar,  mathematically, 
for  most  systems.  Some  judgement  as  to  the  oscillatory  contribution  of  a 
mode  can  then  be  made  once  the  characteristic  damping  and  oscillation 
frequencies  of  the  mode  are  known. 
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